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We describe the analytic calculation of the fermionic two-loop QCD corrections to the heavy-quark pair production 
process in the quark-antiquark channel. 
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<n: 1. introduction 



' The top quark is the heaviest fermion of the Standard Model. Since its discovery at the Fermilab Tevatron 
its mass has been measured to within a few percent, while its production cross section and couplings are currently 
known with larger uncertainty. With the large number of top quarks expected to be produced at the LHC, the study 
of its properties will become precision physics. To interpret these upcoming precision data, equally precise theoretical 
^ i' predictions are mandatory. These demand foremost the calculation of higher order corrections in perturbative QCD. 

• At present, the top quark pair production cross section is known to next-to-leading order (NLO) in the QCD 
O , coupling constant [2]. For this process, the resummation of next-to- leading logarithmically enhanced corrections 
(NLL) improves upon the fixed-order NLO prediction [3]. Electroweak one-loop corrections to ti production are 
equally available For the top quark pair production cross section, which is expected to be measured to within 
a few percent accuracy, the currently available theoretical prediction is not sufficiently precise. Recent studies j3] 
00 indicate a scale uncertainty on these predictions of 7%, and a parton distribution uncertainty of 6%. While the latter 
' may be improved upon by more precise determinations of the parton distribution functions at HERA and LHC, the 
. former requires the calculation of perturbative corrections at next-to-next-to-leading order (NNLO) in QCD. 

The calculation of the full NNLO corrections to the top quark pair production cross section requires three types of 
ingredients: two-loop matrix elements for qq ti and gg — > tt, one-loop matrix elements for hadronic production of 
tt+{l parton) and tree-level matrix elements for hadronic production of tt+(2 partons). The latter two ingredients 
were computed previously in the context of the NLO corrections to tt-|-jet production [0]. They contribute to the tt 
. ^ ' production cross section through configurations where up to two final state partons can be unresolved (coUinear or 
soft), and their implementation thus may require further developments of subtraction techniques at NNLO. 

Both two- loop matrix elements were computed analytically in the small-mass expansion limit s, \t\, \u\ ^ in 
starting from the previously known massless two-loop matrix elements for qij q'q' Q and gg — > qq jpl. An exact 
numerical representation of the two- loop matrix element qq — > ti has been obtained very recently [lO| . In 11 1 
we computed all two-loop contributions to qq — > ti arising from closed fermion loops in a compact analytic form, 
providing also a first independent validation of the results of Our results allow for a fast numerical evaluation 

and permit the analytical study of the cross section near threshold. In the rest of this proceeding we briefly discuss 
the structure of the two-loop fermionic corrections and the calculational techniques employed to evaluate them. 
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2. STRUCTURE 

The scattering process we consider is q{pi) + q(ji2) ^ t{p3) + t{p4) in Euclidean kinematics, where = for z = 1, 2 
and p| = —m^ for i = 3,4. The Mandelstam variables are defined as follows: s = — (pi +^2)^, t = — {pi — Ps)^ , 
u = — (pi — Pi) . Conservation of momentum implies that s + t + u = 2m^. The squared matrix element (averaged 
over the spin and color of the incoming quarks and summed over the spin of the outgoing ones), calculated in 



XXXIV International Conference on High Energy Physics (ICHEP08) 29 Jul - 5 Aug 2008, Philadelphia PA, USA 



d = A — 2e dimensions, can be expanded in powers of the strong coupling constant as as follows: 

2 



\M\ {s,t,m,e) 
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A2 + 0\ 



The tree-level amplitude involves a single diagram and its contribution to Eq. ([1]) is given by 

' (t — + {u — 2m? 



Ao = 47V, 



— e 



(1) 



(2) 



where Nc is the number of colors and Cp = {N^ — l)/2A'c. ^ 
The NLO term Ai in Eq. ^ arises from the interference of one- loop diagrams with the tree-level amplitude 2]. 
The NNLO term A2 consists of two parts, the interference of two-loop diagrams with the Born amplitude and the 
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The latter term ^l^^^-* was studied 



interference of one-loop diagrams among themselves: A2 — ^2 "i" -^2 • -^^^"^ lauuci ucim ^2 
extensively in A^^^'' can be decomposed according to color and flavor structures as follows 
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where Ni and Nf^ are the number of light- and heavy-quark flavors, respectively. The coefficients A, B, . . . , in 
Eq. ([3]) are functions of s, t, m, and e. These quantities were calculated in Q in the approximation s, \t\, \u\ 3> rr?. 

(2x0) 

For a fully differential description of top quark pair production at NNLO, the complete mass dependence of A2 is 
required. An exact numerical expression for it has been obtained in 10|. In 11| , we derived exact analytic expressions 
for all the terms in Eq. ([3]) arising from two-loop diagrams involving at least a fermion loop (i.e. the coefficients 
Di, Ei, Fj with i = l,h and j = h, Ih), providing also an independent confirmation of the results of HQ. 



3. CALCULATION 



The two-loop Feynman diagrams for qq tt were generated with QGRAF 13|. The interference with the tree- 
level amplitude, as well as the color and Dirac algebra, were simplified by using a FORM IJ] code. Out of the ^ 
200 two-loop diagrams contributing to the amplitude, about 60 are proportional to Ni and/or N^- There is only 
one two-loop box topology contributing to the Ni part of the squared amplitude, and a single other two-loop box 
topology proportional to N^. These two box topologies are very similar to the ones encountered in the evaluation of 
the two-loop QED corrections to Bhabha scattering ,[3], and can be evaluated with the same techniques. 

All two-loop integrals appearing in these amplitudes are reduced to a set of master integrals (Mis) using two 
independent implementations of the Laporta algorithm [l3| . Only part of these Mis were available in the literature [l^ 
from previous two-loop calculations of the heavy quark form factors fiot and amplitudes for Bhabha scattering 
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The remaining Mis were evaluated in [ill by employing the differential equation method 



All the Mis were calculated in the non-physical region s < 0. The transcendental functions appearing in the Mis 
are one- and two-dimensional harmonic polylogarithms (HPLs) [22j of maximum weight four and three, respectively. 
Both sets of functions can be rewritten in terms of conventional Nielsen's polylogarithms. 

Following the procedure outlined in the present section, it was possible to obtain the expression of the bare 
squared matrix elements involving diagrams proportional to Ni and/or N^- The UV divergencies were renormalized 



ll|. In order to cross check our analytical results, we expanded them 

the second 



in a mixed scheme described in detail in 

in the s, |u| ^ m? limit. The first term in the expansion agrees with the results published in 



order term agrees with the results found in the Mathematica files included in the arXiv version of [l0|. We also find 
complete agreement with the numerical result of Table 3 in [l^] , corresponding to a phase space point in which the 
s, \t\, \u\ S> m? approximation cannot be applied. 
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